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SUMMARY 


Various  processes  are  discussed  that  affect  the  electrical  coupling  of 
the  ionosphere  and  the  magnetosphere  and  the  formation  and  distribution 
of  plasma  irregularities.  The  properties  of  electrostatic  potential  dif- 
ferences along  magnetic  field  lines  are  inferred  from  analyses  of  the 
pitch-angle  and  energy  distributions  of  the  electrons  and  ions  measured 
with  the  S3-3  satellite.  Such  potential  differences,  as  high  as  8 kV, 
appear  about  70%  of  the  time  along  field  lines  at  latitudes  *65°.  In 
a detailed  analysis  of  an  event  observed  at  an  altitude  of  7280  km,  in- 
variant latitude  -70°,  1800  MLT,  1054  UT  on  15  September  1976,  the  poten- 
tial differences  below  and  above  the  satellite  were  deduced  from  the 
electron-distribution  measurements.  The  total  potential  difference  along 
the  magnetic  field  was  found  to  have  an  inverted-V  type  structure, 
increasing  along  the  satellite  path  from  zero  to  a maximum  of  about  4.3 
kV  and  back  to  zero.  The  spatial  extent  of  the  structure  in  the  north- 
south  direction  is  about  340  km  if  it  is  assumed  that  the  structure  is 
stationary.  Furthermore,  the  potential  distribution  above  the  satellite 
was  estimated  by  computing  the  number  densities  of  the  electrons  and  ions 
from  the  measured  fluxes  and  using  the  quasi-neutrality  condition.  The 
corresponding  electric  field  was  found  to  increase  from  zero  at  a distance 
of  about  1460  km  above  the  satellite  to  about  0.82  mV/m  at  the  satellite. 
These  results  together  with  observations  in  that  region  of  electrostatic 
waves  with  AC  electric  fields  of  -10  mV/m  indicate  that  the  potentials 
are  supported  along  the  magnetic  field  by  anomalous  resistivity. 

Data  from  the  WIDEBAND  satellite  program  were  investigated  with  the  intent 
of  correlating  scintillation  activity  with  electron  precipitation  events 
and  large-scale  electric  fields.  The  data  obtained  at  mid-latitudes  and 
at  high  latitudes  are  discussed.  Certain  types  of  scintillations  are  identi- 
fied, and  the  morphology  of  the  scintillations  in  the  auroral  regions  is 
described.  In  general,  scintillation  activity  has  a spatial  structure  and 
occurrence  similar  to  discrete  auroral  arcs  and  electric  fields.  Tentative 
interpretations  are  given  for  some  of  the  observations. 


It  was  found  that  widespread  plasma  irregularities  may  result  from  a high- 
altitude  nuclear  explosion.  The  electric  fields  and  currents  in  the 
ionosphere-magnetosphere  system  resulting  from  the  injection  of  debris 
plasma  and  fission  beta-decay  electrons  were  computed  using  a previously- 
developed  model  for  the  debris  motion.  The  currents  flow  along  the  magnetic 
field  lines  at  high  altitudes  and  across  magnetic  field  lines  in  the  iono- 
sphere. For  bursts  in  the  megaton  range,  it  was  found  that  both  the  field- 
aligned  and  cross-field  currents  become  unstable  and  produce  plasma  irregu- 
larities. The  pertinent  plasma- instability  criteria  indicate  that  for 
bursts  in  the  northern  hemisphere  the  irregularities  may  occur  over  a wide 
region  north  and  east  of  the  burst.  To  the  east  the  region  may  extend 
over  a longitudinal  sector  approaching  60° ; to  the  north  the  region  may 
reach  the  auroral  zone. 
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Section  1 
INTRODUCTION 


The  degradation  of  radar  and  communications  systems  resulting  from  plasma 
irregularities  in  the  ionosphere  continues  to  pose  a severe  defense  problem. 
Steady  progress  is  being  made  on  the  understanding  of  the  interference  ef- 
fects and  on  the  prediction  of  the  irregularity  patterns  due  to  naturally- 
occurring  ionospheric  disturbances . Several  experimental  programs , includ- 
ing the  WIDEBAND  experiment,  are  providing  valuable  information  on  the 
effects,  and  theoretical  studies  have  identified  several  mechanisms  that 
are  responsible  for  some  of  the  observed  irregularities.  The  more  perplex- 
ing problem,  of  course,  is  the  prediction  of  the  effects  due  to  the  irregu- 
larities resulting  from  a nuclear  exchange.  Some  work  has  been  done  in 
that  area.  The  known  instability  mechanisms  have  been  applied  to  the  prob- 
lem. However,  much  more  work  is  required.  This  report  is  on  topics  that 
are  pertinent  to  the  occurrence  of  irregularities  in  both  the  natural  and 
nuclear- induced  environments . 

The  electrical  coupling  of  the  ionosphere  and  the  magnetosphere  has  an  im- 
portant influence  on  the  formation  and  distribution  of  the  irregularities. 
It  affects  the  distribution  of  charged  particles  in  the  coupled  regions, 
the  energy  deposition  in  the  atmosphere,  and  the  motions  of  the  magnetic 
tubes  that  contain  enhanced  and  striated  plasma.  The  S3-3  satellite  con- 
tinues to  provide  the  best  available  information  on  this  problem.  Our 
analysis  of  these  data  are  described  in  Section  2.  The  data  clearly  indi- 
cate that  for  conditions  above  auroral  arcs  the  electrical  resistivity 
along  the  magnetic  field  is  often  many  orders  of  magnitude  ( 10^)  greater 
than  the  resistivity  inferred  from  binary  collisions.  Hence,  the  usual 
assumption  made  in  magnetospheric  physics,  that  field  lines  are  equipoten- 
tials,  is  not  always  valid.  Potentials  along  the  magnetic  field  as  high 
as  8 kV  are  inferred  from  the  measurements  of  the  charged-particle  fluxes 
obtained  with  S3-3  satellite.  It  appears  that  these  potentials  are  due 
to  plasma  instabilities  driven  by  field-aligned  currents,  but  the  data  are 
not  in  agreement  with  the  available  theories.  Much  more  information  on 
this  phenomenon  is  required  in  order  to  predict  the  effects  which  would 
result  from  a high-altitude  nuclear  explosion. 
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The  following  section,  3,  describes  our  efforts  to  correlate  the  regions 
of  scintillation  activity  with  precipitation  regions  and  large-scale  elec- 
tric fields.  A sufficient  amount  of  data  from  the  WIDEBAND  program  are 
available  to  enable  us  to  construct  a picture  of  the  general  scintillation 
morphology  in  the  auroral  regions.  In  particular,  we  believe  we  can  identify 
certain  distinct  regions  of  scintillation  activity  with  the  regions  of 
auroral  precipitation.  The  spatial  structure  and  occurrence  of  these 
regions  is  similar  to  the  structure  and  occurrence  of  discrete  auroral  arcs 
and  electric  field  regions.  It  was  not  possible  within  the  scope  of  the 
present  work  to  make  any  detailed  correlations,  but  the  evidence  is  suffi- 
cient to  indicate  that  the  scintillation  regions  identified  here  are 
promising  objects  for  further  study. 

The  hot  plasma  - the  nuclear  debris  and  relativistic  electrons  - injected 
in  the  magnetosphere  by  a high-altitude  nuclear  explosion  (H.A.N.E.)  pro- 
duces widespread  electric  fields  and  currents , both  along  and  across  the 
magnetic  field,  in  the  ionosphere-magnetosphere  system.  In  our  previous 
report  (Ref.  1) , the  electric  field  induced  by  the  atmospheric  heave  and 
the  resulting  current  in  the  E-region  of  the  ionosphere  were  estimated. 

The  current  was  found  to  be  unstable  to  the  production  of  electrostatic 
waves,  according  to  the  Farley-Buneman  mechanism.  This  instability  forms 
magnetic-field  aligned  plasma  irregularities;  but,  for  this  case,  the 
instability  is  limited  to  distances  in  the  ionosphere  comparable  to  the 
dimensions  of  "heave"  region.  In  Section  4 the  electric  fields  and  currents 
due  to  charge-separation  processes  in  the  magnetosphere  resulting  from  the 
hot  plasma  injection  are  estimated.  It  is  found  that  these  currents  are 
unstable  in  the  same  way,  but  over  a much  wider  region.  For  bursts  in  the 
megaton  range,  the  field  aligned  irregularities  may  extend  north  of  the 
burst  point  to  latitudes  approaching  the  auroral  regions  and  east  of  the 
burst  to  a longitudinal  extent  of  about  60  deg. 
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Until  recently  it  was  generally  thought  that  quasi-steady  electric  fields 

could  not  be  supported  along  magnetic  field  lines  at  altitudes  above  a 

few  hundred  kilometers.  At  those  altitudes  the  collision  frequencies  of 

the  charged  particles  with  the  neutral  constituents  of  the  atmosphere  are 

very  small.  The  mobilities  of  the  charged  particles  along  the  magnetic 

field  were  therefore  thought  to  be  so  high  that  the  particles  could  rapidly 

become  redistributed  along  the  field  in  such  a way  as  to  neutralize  any 

electric  field  that  would  tend  to  develop.  However,  as  discussed  in  Ref.  1, 

several  measurements  of  charged-particle  fluxes  in  the  auroral  zones  and 

a few  observations  of  shaped-charge  barium  releases  indicated  the  presence 

of  high  electric  potentials  along  the  magnetic  field.  At  the  present 

time  the  S3-3  satellite  is  providing  valuable  information  on  the  nature 

of  these  potentials.  They  are  observed  in  about  half  of  the  satellite 

passes  at  high  latitudes,  generally  above  auroral  arcs  (Ref.  2).  Quite 

often  these  potentials  are  accompanied  by  field-aligned  currents  greater 
“6  2 

than  10  amps/m  , broadband  AC  electric  fields  of  tens  of  mV/m  associated 
with  electrostatic  ion  cyclotron  waves  (Ref.  3),  and  ions  of  ionospheric 
origin  with  energies  of  several  keV  (Ref.  4).  These  conditions  indicate 
an  anomalous  resistivity  along  the  magnetic  field  of  about  100  ohm-m 

3 

(plasma  density  ~100/cm  ) compared  with  the  classical  resistivity  of  about 
1.5  x 10  ^ ohm-m  (Ref.  5). 


Falthammar  has  reviewed  recently  the  theories  that  have  been  proposed  to 
explain  the  field-aligned  potentials  (Ref.  6).  Two  of  these  theories  pre- 
dict conditions  which  are  in  agreement  with  the  data.  Kindel  and  Kennel 

(Ref.  7)  have  shown  that  a field-aligned  current  with  an  intensity  greater 
“6  2 

than  about  10  Amps/m  and  with  a drift  velocity  that  exceeds  a certain 
minimum  value  is  unstable  to  the  production  of  electrostatic  waves.  More- 
over, they  showed  that  the  instability  favors  the  production  of  electro- 
static ion  cyclotron  waves  rather  than  acoustic  waves;  and,  indeed,  the 
data  reveal  that  the  electrostatic  waves  consist  principally  of  the  ion- 
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cyclotron  waves.  The  interaction  of  the  electrons  carrying  the  current 
with  the  AC  electric  fields  of  the  waves  retards  the  electrons  in  a manner 
similar  to  collisions  in  a resistor.  Accordingly,  a potential  builds  up 
across  this  instability,  or  "anomalous  resistivity",  region.  The  DC  elec- 
tric field  along  the  magnetic  field  must  be  smaller  than  the  turbulent 
AC  wave  field  (Ref.  6).  Hence,  since  the  observed  wave  field  is  typi- 
cally 1-10  mV/m  (Ref.  8),  the  expected  DC  electric  field  is  of  the  order 
of  1 mV/m.  The  electric  field  due  to  anomalous  resistivity  must  therefore 
extend  several  thousand  kilometers  along  the  magnetic  field  in  order  to 
account  for  the  inferred  potentials. 

The  other  theory  which  is  sometimes  supported  by  the  data  was  proposed  by 
Block  and  Falthammar  (Ref.  9 ) . They  suggested  that  an  electric  double 
layer  would  result  from  a field-aligned  current  instability.  The  spatial 
extent  of  a double  layer  is  small,  being  scaled  in  Debye  lengths,  but  the 
electric  field  in  the  double  layer  is  high.  The  measurements  of  Mozer 
et  al,  (Ref.  3)  indicate  the  presence  of  high  DC  electric  field  intensities 
(,v.5V/m)  both  along  and  across  magnetic  field  lines,  but  the  distribution 
of  the  field  resembles  the  field  expected  for  the  paired  oblique  electro- 
static-shock structure  that  has  been  proposed  by  Swift  (Ref.  10).  This 
is  a modification  of  the  double  layer.  The  paired  shocks  are  oblique  to 
the  magnetic  field,  and  the  electric  field  in  the  shock  is  more  extensive, 
being  scaled  in  ion  gyroradii  rather  than  Debye  lengths. 

In  Section  2.2  the  electron  and  ion  fluxes  observed  with  the  S3-3  satellite 
through  one  of  the  field-aligned  potential  regions  are  discussed,  and  the 
analysis  of  the  data  that  yields,  the  potential  differences  below  and  above 
the  satellite  is  described.  In  Section  2.3  the  particle  distributions 
based  on  the  observations  are  transformed  to  altitudes  above  the  spacecraft, 
assuming  the  particle  motions  to  be  adiabatic  in  the  electric  and  magnetic 
fields,  and  the  number  densities  of  the  particles  are  computed  for  assumed 
electric  potential  distributions.  The  potential  distribution  given  by  the 
quasi-neutrality  condition  is  then  determined  and  compared  with  the  distri- 
butions predicted  by  theory. 
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2.2  ELECTRIC  POTENTIALS  INFERRED  FROM  ELECTRON  MEASUREMENTS 

In  Ref.  1 we  described  the  electrostatic  potential  differences  along 
magnetic  field  lines  inferred  from  an  analysis  of  the  pitch-angle  and 
energy  distributions  of  the  electrons  observed  with  the  S3-3  satellite. 

The  analysis  was  for  an  event  observed  at  1054  UT  on  15  September  1976 
at  18  hr  MLT,  70°  invariant  latitude,  and  at  an  altitude  of  7280  km. 

In  reference  1 the  analysis  was  done  with  an  incomplete  set  of  the  data 
on  the  event.  The  further  data  has  since  been  obtained.  The  results 
described  below  are  based  on  an  improved  analysis  of  all  of  the  data 
available  on  this  event.  Moreover,  the  analysis  yields  the  distance 
over  which  the  electrostatic  potential  extends  above  the  spacecraft. 

The  motion  of  the  satellite  through  the  region  sustaining  the  electric 
field  is  depicted  in  Figure  1.  The  satellite  moved  through  this  region 
toward  higher  latitudes.  On  successive  spin  periods,  the  satellite  was 
located  at  the  positions  indicated  by  the  numbers  1-4  in  the  figure. 

The  energy  ranges  of  the  four  detectors  used  to  observe  the  electrons 
are  listed  in  Table  1. 

Table  1 

Energy  Ranges  of  Electron  Detectors 

Detector 

CME  A 
CME  B 
CME  C 
CME  D 

The  counting  rates  of  the  detectors,  on  the  various  spin  periods,  as  func- 
tions of  the  pitch-angle  of  the  electrons  are  shown  in  Figures  2-5.  Note 
that  in  Figures  2-4,  the  ordinates  have  been  displaced  by  a decade  to  sepa- 
rate the  curves.  The  pitch-angle  notation  is  chosen  such  that  electrons 
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Figure 


1 Schematic  diagram  of  motion  of  satellite  through 

region  containing  electrostatic  potential  (see  text). 
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having  pitch  angles  (a)  of  0 and  180°  move  along  the  magnetic  field 
toward  and  away  from  the  near  atmosphere,  respectively.  The  numbers 
at  the  curves  denote  the  spin  number,  and  the  superscripts  - and  + denote 
the  1st  and  2nd  pitch-angle  scans  that  are  made  during  each  spin  period. 
Since  the  counting  rates  of  the  highest-energy  detector,  CME  D,  were  low, 
the  counting  rates  shown  in  Figure  5 were  averaged  over  the  separate 
scans  of  the  spin  periods. 

On  spin  1 note  that  the  counting  rates  of  all  of  the  detectors  indicate 
that  the  electron  flux  was  very  nearly  isotropic  except  for  the  loss  cone 
due  to  the  near  atmosphere.  The  counting  rates  of  the  detectors  had  been 
the  same  for  many  spin  periods  preceding  spin  1,  and  the  electron  flux 
was  very  closely  represented  by  a Maxwellian  distribution  with  a number 
density  of  0.76  and  a temperature  of  0.8  keV. 


On  the  succeeding  spins  the  presence  of  the  potential  below  and  above  the 

satellite  is  more  clearly  demonstrated  by  the  CME  B and  CME  C counting 

rates.  In  Figures  3 and  4 note  that  the  loss  cones  become  progressively 

wider  and  deeper  as  the  satellite  moves  from  its  position  on  spin  1 to 

that  on  spin  3.  The  edge  of  the  loss  cone,  a , is  related  to  the  poten- 

c 

tial  beneath  the  satellite  through  the  conservation  of  the  magnetic  moment 
of  the  electrons,  i.e., 


or 


w sin  a [w  -+-  q ( d)  - <t>  ) 1 

s c s n m s 


-1,B  . (%  " 9s),  ,1/2 

+ q- *-)  ] 


a = sin  [—  (1  + q- 
c B w 

t s 


(1) 

(2) 


Here,  wg  and  B^  are  the  electron  energy  and  the  magnetic  field  intensity 

at  the  satellite,  q is  the  absolute  value  of  the  electron  charge,  and 

<t>m  -<i>s  is  potential  difference  between  the  satellite  and  limiting  point 

B = Bt  where  the  electron  mirrors  above  the  atmosphere  before  it  suffers 

substantial  collisions.  Note  in  particular  from  Eq.  (2)  that  for  a fixed 

potential  difference  the  loss-cone  width,  w-  a , increases  as  w becomes 

c s 

smaller.  This  effect  is  clearly  demonstrated  by  the  data  in  Figures  3 and  4. 
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31 


PITCH  ANGLE  (DEG.) 


Figure  2 


Counting  rate  of  CME  A detector  as  function  of 
pitch  angle.  The  numbers  at  the  curves  denote 
the  satellite  spin  number,  and  the  superscripts 

- and  + denote  the  preceding  and  succeeding  pitch-angle 
scans  made  during  the  satellite  rotation. 
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Figure  3 Counting  rate  of  CME  B detector  as  function  of  pitch 
angle.  (See  caption  of  Figure  1.) 
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Figure  5 Pitch-angle  data  obtained  on  spins  1 and  3 by  the  CME  D detector. 


The  loss  cones  become  deeper  as  <p  - <p  increases  because  the  backscattered 

m s 

electrons  have  lower  energies  than  the  primaries,  hence  fewer  of  them  can 
penetrate  the  retarding  potential  in  order  to  reach  the  satellite. 

The  "butterfly"  pitch-angle  distributions  on  spin  4 indicate  the  presence 
of  a potential  difference  above  the  satellite.  The  analysis  described 
below  indicates  that  as  the  spacecraft  began  its  fourth  spin  (at  the  pitch 
angle  a = 0)  the  potential  above  the  spacecraft,  <j>s , had  increased  to  about 
0.8  kV  and  the  field  ratio  B^/Bg  was  about  0.77.  This  potential  difference, 
by  accelerating  the  electrons  along  the  field,  caused  the  pitch  angles  of 
the  electrons  to  decrease,  and  the  effect  was  more  pronounced  for  the 
lower  energy  electrons.  Accordingly,  as  shown  in  Figure  4,  the  counting 
rate  of  CME  C was  relatively  constant  from  a = 0°  to  60°  on  scan  4 . All 
of  the  electrons  within  the  energy  range  of  the  detector  were  detected  at 
those  pitch  angles.  As  the  detector  continued  to  rotate  toward  a = 90°, 
it  detected  progressively  lower  fluxes  because  only  the  higher-energy 
electrons  of  the  spectrum  had  the  higher  pitch  angles.  As  the  detector 
rotated  from  a = 90°  to  a =140°,  it  observed  principally  the  electrons 
which  were  reflected  by  the  magnetic  field  below  the  spacecraft;  these 
electrons  were  primaries  but  with  an  augmented  lower-energy  component  due 
to  trapping  between  the  electric  field  above  the  spacecraft  and  the  mag- 
netic field  below.  On  the  return  scan  (4+)  the  detector  observed  a higher 
fraction  of  the  trapped  electrons  among  the  down-going  electrons  (a  < 90°). 

The  counting  rate  of  the  CME  B detector  on  spin  4 is  due  to  different  frac- 
tions of  the  electron  sources.  This  detector  responded  to  electrons  in 
the  range  0.352  to  1.13  keV.  Hence,  the  potential  above  the  spacecraft 
accelerated  most  of  the  primary  electrons  to  energies  which  could  not  be 
observed  by  the  detector.  Here,  the  peaks  in  the  pitch-angle  distribution 
are  due  to  the  higher  energy  electrons  within  the  sensitivity  of  the  detec- 
tor which  are  principally  degraded  primaries  due  to  scattering  in  the 
atmosphere.  Toward  the  minimum  at  a = 90°,  the  counting  rates  are  due  to 
progressively  lower-energy  electrons  of  the  atmospheric  backscatter.  Again, 
the  "symmetry"  about  90°  is  due  to  the  reflection  of  the  electrons  by  the 
electric  field  above  the  spacecraft.  A high  fraction  of  these  lower  energy 

electrons  are  trapped  locally. 
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The  counting  rate  of  detector  CME  A (See  Fig.  2)  beyond  the  second  half  of  scan 
2 and  just  prior  to  the  end  of  scan  4+  is  principally  due  to  locally-trapped 
and  untrapped  backscattered  electrons.  During  that  time  the  potential  difference 
above  the  satellite  exceeded  the  high-energy  limit  of  the  detector  response 
(.238  keV) . Hence,  the  primary  electrons  were  accelerated  to  energies  beyond 
the  sensitivity  region  of  the  detector. 


The  detector  CME  D is  not  too  useful  in  the  analysis  of  this  event  because 
its  response  is  to  energies  appreciably  above  the  magnitudes  of  the  poten- 
tials . 

The  energy  ranges  of  the  detectors  (see  Table  1)  are  too  broad  to  resolve 

in  sufficient  detail  the  changes  in  the  spectrum  due  to  the  potential 

and  pitch-angle  changes.  Therefore,  the  following  procedure  was  used  to 

determine  the  potentials.  The  electron  flux  above  the  satellite,  where 

the  potential,  $ , was  zero  and  B=B^  (see  Fig.  l),was  assumed  to  be  the 

same  as  that  measured  on  spin  1 at  B=Bg . This  assumption  seems  plausible 

because,  as  mentioned  above,  the  electron  flux  measured  on  spin  1 was  the 

same  as  it  had  been  for  many  of  the  previous  spin  periods.  A further 

assumption  was  that  the  electrostatic  potential  increased  monotonically 

from  | = 0 at  B.  to  tp  = 0 below  the  satellite,  and  remained  constant, 

1 m 

at  $ , to  the  top  of  the  atmosphere  (300  km)  where  B=B  (see  Fig.  1).  The 
m t 

flux  at  B=B^  was  then  transformed  to  the  position  of  the  satellite  using 

Liouville's  theorem  for  various  values  of  B,  , <t>  , and  4>  . This  transformed 

Is  m 

flux  included  the  downgoing  electrons  which  arrived  at  the  satellite  directly 
from  B^ , the  upgoing  electrons  which  had  mirrored  below  the  satellite  but 
above  the  "top"  of  the  atmosphere,  and  the  downgoing  electrons  which  were 
reflected  above  the  satellite  by  the  repulsive  potential  - <j)g. 

An  iteration  procedure  was  used  to  compute  the  equilibrium  flux  at  the 

satellite  due  to  the  atmospheric  backscatter  for  various  values  of  B^ , <f>g, 

and  (f  . First,  the  flux  at  B.  was  transformed  to  B where  it  was  used  as 
m It 

input  to  the  AURORA  code  to  compute  the  backscattered  flux.  Those  electrons 
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which  were  reflected  by  the  potential  - $ were  then  added  to  the 

m 

previously  determined  flux  of  electrons  incident  on  the  atmosphere,  and 
the  AURORA  code  was  used  again  to  compute  the  enhanced  backs cattered  flux. 
The  computations  were  continued  in  this  manner  until  the  backscattered 
flux  converged  to  a constant  value.  Owing  to  the  repulsive  potential,  -<j>  , 
above  the  satellite,  the  backscattered  electrons  contributed  to  the  down- 
going fluxes  at  the  satellite  as  well  as  to  the  upgoing  fluxes. 

The  theoretical  fluxes  at  the  satellite  were  then  integrated  over  the 
energy-response  functions  of  the  detectors,  for  various  values  of  the  elec- 
tron pitch  angles,  and  the  resulting  counting  rates  were  compared  with  the 
experimental  values.  Those  values  of  the  electrostatic  potentials  below 
and  above  the  spacecraft  which  brought  the  theoretical  counting  rates  into 
agreement  with  the  experimental  values  are  shown  in  Figure  6.  On  spin  5 
the  electron  and  energetic- ion  fluxes  dropped  off  suddenly  to  the  noise 
levels  of  the  detectors  indicating  an  absence  of  the  potential.  Hence, 
the  analysis  reveals  an  inverted-V  type  structure  similar  to  the  paired, 
oblique  electrostatic  shock  structure  proposed  by  Swift  (Ref. 10  ) to 
explain  various  phenomena  which  have  been  observed  in  and  above  the  auroral 
zones.  However,  important  differences  exist.  The  structure  shown  in 
Fig.  6,  when  transformed  along  field  lines  to  an  altitude  of  100  km, 
has  a north-south  extent  of  about  70  km.  This  is  much  wider  than  structure 
proposed  by  Swift,  which  is  comparable  to  the  width  of  a discrete  auroral 
arc.  The  observed  width  is  closer  to  that  inferred  by  Frank  and  Ackerson 
(Ref.  11)  to  explain  the  inverted-V  shaped  electron  fluxes  that  they  ob- 
served with  the  Injun  5 satellite  at  high  latitudes.  The  energies  of  the 
electron  fluxes  observed  by  these  experimenters  at  the  apex  of  the  struc- 
ture are  also  consistent  with  the  peak  value  of  the  potential  shown  in 
Fig.  6.  Furthermore,  the  electric  field  along  the  magnetic  field  appears 
to  be  much  smaller  than  that  required  by  Swift's  distribution:  the  mean 
value  of  the  electric  field  can  be  estimated  from  the  Spin  4 data  which 
indicate  a potential  above  the  spacecraft  of  0.8  kV  over  a magnetic-field 
ratio  B^/B^  of  0.77  (distance  of  1460  km).  These  values  imply  a mean 
field  of  0.8  kV/1460  km  or  ~0.5  mV/m.  The  analysis  described  in  Section 
2.3  provides  further  information  about  this  electric  field. 
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Figure  6 
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Electrostatic  potentials  along  magnetic  field  inferred  from  pitch-angle 
and  energy  distributions  of  electrons  measured  on  spins  1 through  4.  The 
circles  denote  the  best  estimates  of  the  potential  difference  below  the 
satellite.  The  uncertainties  of  these  values  are  indicated  by  the  error 
bars.  The  triangles  denote  the  potential  difference  below  the  satellite 
inferred  from  the  depths  of  the  loss  cones  alone.  The  squares  denote 
the  total  potential  difference,  which  includes  the  potential  difference 
above  the  satellite. 
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As  mentioned  in  the  introduction,  in  some  of  these  events  the  direct 
electric-field  measurements  indicate  the  presence  of  potential  structures 
that  are  more  nearly  comparable  to  those  described  by  Swift  and  with 
electric  fields  of  hundreds  of  mV/m  both  along  and  across  magnetic  field 
lines  (Ref.  3). 

On  spins  2,  3,  and  4 high  fluxes  of  H+  and  0+  ions  moving  upward,  closely 
aligned  with  the  magnetic  field  (a  ^170°),  were  also  observed.  The 
acceleration  of  ions  from  the  ionosphere  by  the  inferred  potentials  is 
expected.  However,  on  this  event  the  ion  data  do  not  yield  additional 
information  on  the  magnitudes  of  the  potentials.  Several  complications 
can  be  mentioned.  (i)  The  ions  evidently  also  interact  strongly  with 
the  electrostatic  ion  cyclotron  waves  which  seem  always  to  be  present 
during  these  events  (Ref.  13).  (ii)  The  travel  time  of  the  ions  along 
the  acceleration  region  is  appreciable  (*  10  sec  for  a 1 keV  proton) , 
hence  the  ion  spectrum  does  not  reveal  the  "instantaneous"  conditions 
below  the  spacecraft.  (According  to  Figure  6,  the  potentials  change 
appreciably  along  the  satellite  path  in  10  sec.)  And  (iii)  the  ion 
spectra  are  not  well  resolved  in  time;  each  of  the  three  ion  spectrom- 
eters samples  a specific  energy  step  for  16  sec. 

On  this  event,  the  0+  ion  fluxes  were  much  lower  than  the  H+  fluxes. 

The  spectrum  of  the  observed  proton  flux,  the  directional  flux  averaged 
over  the  "source  cone"  (170°  s a i 180°)  is  shown  in  Fig.  7.  The 
numbers  at  the  data  points  are  the  spin  periods  over  which  the  data  were 
acquired.  The  open  triangles  are  the  data  that  were  obtained  on  spin  4 
by  researchers  of  the  Aerospace  Corporation  (Ref.  14)  who  had  an  ion 
detector  on  the  same  satellite.  The  curves  drawn  in  the  figure  are  the 
fluxes  corresponding  to  Maxwellian  distributions  but  only  within  the 
very  restrictive  pitch-angle  range;  the  solid  curve  is  for  a "temperature 
of  0.2  keV  and  the  broken  curves  ia  for  a "temperature "of  0.4  keV.  TTie 
data  seem  to  indicate  that  the  protons  are  suffering  anomalous  collisions 
while  being  accelerated  along  the  field,  and  perhaps  also  that  "runaway" 
protons  are  being  accelerated  throughout  the  potential  region. 


ENERGY  (KeVl 

Directional  flux  of  protons  averaged  over  the 
10-deg  source  cone  as  a function  of  the  proton 
energy.  The  numbers  at  the  points  denote  the 
satellite  spin  periods  during  the  measurements. 

The  triangles  denote  the  data  obtained  by  the 
Aerospace  detector  on  the  same  satellite  (Ref.  14) 


2.3  DISTRIBUTION  OF  ELECTRIC  POTENTIAL  ALONG  MAGNETIC  FIELD 


As  discussed  in  the  introduction,  the  distribution  of  the  potential  along 
the  magnetic  field  provides  a clue  to  the  mechanism  responsible  for  the 
potential.  An  attempt  was  made  to  estimate  the  potential  distribution 
above  the  spacecraft,  on  spin  4,  by  computing  the  number  densities  of  the 
ions  and  electrons  in  the  potential  region. 

The  potential  distribution  is  given  by  a solution  of  Poisson's  equation, 


d2<fr 

ds2 


_ i_ 
k 

o 


<ni 


- n ) 
e 
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where  s is  a distance  along  the  magnetic  field,  q is  the  absolute  value 
of  the  electron  charge,  n^  and  n^  are  the  number  densities  of  the  ions 
and  electrons,  and  is  permittivity  of  free  space.  However,  if  the 
scale  of  the  potential  distribution  is  large,  as  indicated  by  the  analysis 
described  in  the  previous  section,  the  left-hand  side  of  (3)  is  small, 
and  the  potential  distribution  can  be  estimated  by  finding  that  distribu- 
tion that  makes  n.  equal  to  n (see,  e.g.,  Persson,  Ref.  15). 
l e 

The  number  densities  above  the  satellite  (B  < B ) were  computed  for  the 

s 

following  particle  fluxes: 

(I)  the  flux  of  electrons  moving  downward  from  B^  < B which  was 
assumed  to  be  the  same  as  that  observed  on  spin  1, 

(II)  the  observed  flux  of  electrons  moving  upward,  and 

(III)  the  observed  flux  of  protons  moving  upward. 

The  0+  ion  flux  was  neglected  because,  as  mentioned  previously,  it  was 
small  in  comparison  to  the  proton  flux.  Also,  the  magnetospheric  protons 
were  not  included  in  this  analysis;  essentially  no  magnetospheric  protons 
were  observed  within  the  energy  range  of  the  ion  spectrometers. 
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The  general  procedure  was  to  transform  the  fluxes  from  B or  B,  to  the 

s 1 

region  B > B > B by  using  Liouville's  equation, 

S 1 


j (w , a , B)  _ jb(wb*  V V 
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where 
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(5) 


Here,  j is  the  directional  flux  (number/cm  -sec*sr .keV) , the  subscript  b=s 

or  1 refers  to  the  quantities  at  the  satellite  (B=B  ) or  at  B=B,  where  <J>=0, 

s 1 

and  \j>  is  the  energy  gain  or  loss  of  the  particles  that  move  through  the 
potential  from  B^  to  B.  The  pitch  angles  are  related  through  the  conserva- 
tion of  the  magnetic  moment. 
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The  number  density  was  then  computed  by  dividing  J by  the  velocity  v , 

3 

which  yields  the  number  of  particles/cm  • sr  • keV,  and  integrating  over 

1/2 

the  appropriate  ranges  of  a and  w.  Hence,  by  putting  v = (2w/m)  , the 

number  density  is  given  by  the  equation 
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The  number  densities  were  computed  for  a trial  potential  distribution  4>(B), 

and  the  results  for  the  ions  and  the  electrons  were  compared  over  the  full 

range  Bg  >_  B >_  B^.  The  boundaries  of  <j>(B)  were  maintained  at  the  values 

given  by  the  electron  pitch-angle  distributions,  viz.,  <}> ( B^ ) = 0 and 

$ (Bg)  = 0.8  kV.  The  computations  were  then  repeated,  using  a different 

shape  for  4>(B)  until  the  number  densities  of  the  ions  and  electrons  were 
» 

approximately  equal.  The  computations  were  facilitated,  as  discussed  below, 
by  fitting  the  data  with  simple  mathematical  functions. 
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2.3.1  Number  Density  of  Magnetospheric  Electrons 


The  flux  of  electrons  moving  downward  from  is  referred  to  here 
as  the  magnetospheric  electrons  (flux  (I)  listed  above) . As  mentioned  pre- 
viously, this  flux  was  assumed  to  be  the  same  as  that  observed  on  spin  1. 

It  was  very  closely  represented  by  a Maxwellian  velocity  distribution  with 
e 3 

a number  density  n of  0.76/cm  and  a temperature  w^  of  0.8  keV.  The  flux 
j at  B , corresponding  to  the  Maxwellian  distribution,  may  be  written  as 


jl(VVBl)  = nQe  k^w^e  Wl^We,  0 <_  y 

3 3 -1/2 

where  k = (2tt  raw  ) 
e e 

From  Eqs.  (4),  (5),  and  (9),  the  flux  at  B transforms  to 

j (w,a,B)  = n ek  we'(w-*)/we 
J o e 


where 


w = w^  + (B) 


iJj(B)  = q<p (B)  is  the  energy  gained  by  the  electrons.  For  electrons  moving 
downward,  0 ± o.  <_  ir/2,  the  ranges  of  a and  w for  the  integration  (8) 
are  obtained  from  Eq.  (7).  The  are, 


0<a<a  =Sin 
“ ~ H 


_1  If^-w5 


where 


0 < a < tt/2 


= ip/  (1-Bj^/B) 


Hence,  the  number  density  of  the  downgoing  electrons  of  the  flux  I is, 
e i|>/w^  w,  a„  (w) 


t n e 

, 1 , _ o e 

(ne  ’ 3 1/2 

Downgoing  (itw^  ) 


, f -w/w 
+ dw  e e 


1 ■’0 


AT  Sin  ada 
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Downgoing 


e 

n 

' 1-  Erf 

[<*/»  >1/21 

— 

o 

L e J 

2 


+ 2_  (B/B1 

nF 


B/B  _1 
1 


(15) 


Here,  D(x)  is  Dawson's  integral: 

rx 


D(x)  = e 


2 

efc  dt 


(16) 


The  number  density  n1  must  include  those  electrons  of  flux  I that  are 
e 

reflected  at  B",  where  B £ < B and  \p  <_  \p'(B")  < \p  • The  electrons  that 

s s 

are  reflected  at  B > are  included  in  the  flux  (II)  listed  above.  The 
pitch  angle  at  B of  an  electron  that  mirrors  at  B'  is  given  by  the 
equation, 

sin2,  -f.  (1  + ) (I7) 

m B w 7 

The  upper  limit  of  the  pitch-angle  integral  in  (8)  for  the  reflected 
electrons  is  it-  a*(w)  where  a^(w)  is  the  minimum  value  of  a for  a given 
value  of  w.  By  differentiating  (17)  with  respect  to  B'  and  setting  the 
result  equal  to  zero,  the  mirror  point,  B',  corresponding  to  a*,  is  found 
to  be  at, 

(df^B;  = (w  + <T(b;)  - <i>  (b))/b;  as) 

(In  Appendix  2A  at  the  end  of  this  section,  it  is  shown  that  the  extreme 

2 2 

value  of  (17)  is  a minimum  only  if  d <p'/dB'  > 0;  but  only  such  potential 

functions  (positive  curvature)  have  been  found  to  give  appropriate  solu- 
tions.) Mathematical  functions  were  used  for  tJj(B)  , and  a modification 
of  Newton's  method  was  used  to  solve  (18)  for  B'.  The  value  of  a*  was 
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then  computed  by  putting  in  Eq.  (17).  The  minimum  energy,  w^,  given 
by  Eq.  (19)  for  electrons  mirroring  at  B,  is 


= B cfy/dB 


(19) 


Hence,  the  number  density  of  the  upgoing  electrons  of  flux  I is, 


(neI} 


upgoing 


e \h/w 
n e e 
= __o 

/ _ 3.1/2 


dw 


5-0* (w) 

— w / w 

e e r~  sinada 
/w 


(20) 


w*  "/2 

This  contribution  to  the  number  density  was  integrated  numerically.  Its 
magnitude  given  by  the  final  iteration  was  also  verified  analytically  by 
using  a function  tj> (B)  that  closely  approximated  the  final  distribu- 
tion in  the  range  0.8  £ B/Bg  £l  (see  Section  2.3.4),  but  that  also 
allowed  (20)  to  be  integrated  in  closed  form.  This  function  is  given 
by  the  equation, 

4>  = c (B/Bs-a)  (21) 


Here  c = 3.72  kV  and  a *»  0.785.  For  this  potential  distribution, 
5*' 


B.,  = B , w = cB/B  , and  the  integration  of  (20)  gives, 


T nn  r 1/2 

(“.  ) . ~~  f (1  - f > 

upgoing  s 


exp  - 


L"  "e(1-B/0j 


\p  B/B  - \p 
s s 


(22) 


The  final  value  of  n is  the  sum  of  the  components  given  in  (15)  and  (20) 

e i 

or  (22).  Note  that  from  (15)  and  (22),  the  value  of  n at  B,  (ii/=0)  is, 

e 1 


■>i>  - -i 


ip  /w 

s e 

B 1/2  ~ B /B  -1 

i + (i  - r)  e 

s 


(23) 


For  B, /B  = 0.77  and  \p  =0.8  keV,  the  contribution  to  the  number  density 
Is  s 

by  the  reflected  electrons  from  (B^  B B^)  is  less  than  2%.  It  is  small 
because  only  the  high-energy  (w  > 3.72  keV)  electrons  of  the  primary 
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distribution  become  reflected  between  B,  and  B . Note  also  that,  as 

1 s 

expected,  if  the  potential  tp^  were  zero  and  the  magnetic  mirror  were 

very  strong  B^/B  -*■  0,  nearly  all  of  the  downgoing  primaries  would 

^ I e 

be  reflected  and  n^  would  approximately  equal  n^  . 

Toward  higher  B values  n ^ decreases  monotonically . 

e 

2.3.2  Number  Density  of  Observed  Electrons  Moving  Upward 

The  electron  flux  at  the  satellite  with  pitch  angles  in  the  range 
90°  £ a <^180°  was  estimated  by  fitting  the  functions, 


e -ws/Wi  T H 

j.(w,a,B)=n.kwe  , a ? c a < a . 

si  sss  les  si  — s — si 


to  the  data  (Figures  2-5)  in  certain  pitch-angle  intervals.  The  "best 

0 Li  H 

fit"  parameters  n.  and  w.  and  the  pitch-angle  limits  a . and  a . are 
ii  si  si 

listed  in  Table  2. 


Table  2 

PARAMETERS  OF  ELECTRON  FLUX  OBTAINED  BY 
FITTING  SPECTRAL  FORM  GIVEN  BY  EQ.  (23)  TO  DATA 


INDEX 

(cm  ) 

w. 

l 

(keV) 

2-i 

a . 
si 

(Deg.) 

n 

asi 

(Deg.) 

1 

1.55 

0.35 

90 

110 

2 

1.85 

0.45 

110 

120 

3 

2.10 

0.59 

120 

138 

4 

1.43 

0.61 

138 

148 

5 

0.50 

0.70 

148 

158 

6 

0.20 

0.50 

158 

180 

(24) 
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Above  the  spacecraft  the  flux  transforms  to 


/ e -(w  + v) /w.  L H 

j.(w,a,B)  = n.k  we  A 1,  a.  < a < a. 

1 le  i — — 1 


where  w = w - x,  X=  q (4>„“4>) 

o S 


(25) 

(26) 


For  the  upgoing  electrons,  the  pitch-angle  and  energy  limits  given 
by  the  magnetic-moment  conservation  equation,  are 


and 


aL  • _1 

a ^ = sin 


H . -1 
= sm 


2 L 


f (1  + - ) sin  a r 
B w si 

s 


. 2 H 


B (1  + w}  Sin  “si 
s 


1/2 


» w ^.W 


1/2 


, w > w2 


(27) 


(28) 


where , 


X8 


L,H 


w 


1,2  L,H 

1_8i 


(29) 


L,H  B . 2 L,H 

8i  = I sln  asi 
s 


(30) 


For  w < w <_  w , the  pitch-angle  integration  in  (8)  is  from  ir/2toa.. 

1 2 L H ^ 

For  w > w^  it  is  from  to  a^.  By  substituting  (25)  into  (8)  and 

performing  the  integration  over  the  pitch  angle,  the  number  density 

may  be  written  as, 

i 


6 e "X/Wi 

n _ c ni  e 

*e  , 3,1/2 

i=l  (ttw±  ) 


[1-(1+  X JgL  ]1/2  e 1 Jw  dw 

W 1 


4l-(l  + J )gi]1/2  -2[1-(1  + J)gJ  11/2 j e 1/w  dw  + R 


(31) 


where  R^  represents  the  corresponding  integrals  for  the  electrons  that 
are  reflected  by  the  electric  field.  For  w < ip , the  electrons  are 
reflected  regardless  of  their  pitch  angle.  Hence,  for  w*  and  w^  less 
than  \p , R will  include  the  integrals  shown  in  (31)  except  that  the  upper 
limit  of  the  second  integral  will  be  ip . If  only  w*  is  less  than  ip  , R^ 
will  include  only  the  first  integral  in  (31)  with  the  upper  limit  being 
ip  instead  of  w^. 


For  w > ip , the  minimum  pitch  angle ,a^,  of  the  reflected  electrons  will 
be  determined  as  discussed  in  the  previous  section.  If  a is  less  than 


H H L 

n - a.,  the  pitch-angle  integration  for  R.  is  from  n-a..  to  ir-a  . 

f jj  I.  i L 

is  between  ir-a^  and  ir-a  , the  integration  is  from  a^to  n-oi^. 


If  a . 


The  integrals  appearing  in  (31)  can  be  integrated  as  indicated  in  the 
previous  section.  Only  the  part  of  R.  for  w > ip  must  be  evaluated 
numerically.  That  part  of  R^  can  also  be  integrated  by  using  the 
approximate  distribution  given  by  Eq.  (21).  By  using  that  function  for 
$ , the  integration  of  Eq.  (31)  yields, 
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L k WT1-W11  x/w. 

« y K(  --w  ■ — ) Exp( i) 
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- M(t”.  /w.)  Exp(  - -2— i )J  | (32) 

1-aS,, 


where 


L,H 


* aS. 


= if  + 


^ ; S*  = sin-  a 


1-S 


L,H 


L,H 


2 «,L»H 
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(33a) 


K(x) 

M(x) 


1 4 Erf(/x) 


/x  - D(/ x) 


/x  Exp(-x) 


(33b) 

(33c) 
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(33d) 


„L>H 


aS 


= 1L  ( 


L,H 


1-aS1 

L,H 


and,  again,  D(x)  is  Dawson's  integral  (Eq.  16).  This  equation  was  used 

to  verify  the  numerical  values  in  the  region  0.8^.  B/B  £.  1. 

s 


2.3.3  Number  Density  of  Observed  Protons  Moving  Upward 

The  measurements  of  the  proton  flux  were  discussed  in  Section  2.2 
and  the  data  are  shown  in  Figure  7.  The  curves  in  the  figure  are 
given  by  the  function, 

p -w  /w 

j (w  ,a  ,B  ) = n k w e S P,  170°  < a < 180°  (34) 

s s s s ops  — S — 


where 


k 

P 


(2tt~ 


3. 

m w ) 
P P 


-1/2 


(35) 


The  solid-line  curve  is  for  w =0.2  keV  and  the  broken-line  curve 

P 

is  for  w =0.4  keV.  Both  of  the  curves  have  been  normalized  such 
P 

that  number  density  of  the  protons  is  equal  to  the  sum  of  the  number 

densities  of  the  electrons  (n^  + n^)  evaluated  at  B=B  . The  pre- 

6 e s 

dominance  of  the  data,  especially  when  the  Aerospace  data  on  spin  4 
are  included,  indicate  that  the  solid-line  curve  is  more  representa- 
tive of  the  average  flux  in  the  source  cone. 


Above  the  spacecraft  this  flux  transforms  to, 


p 

j (w,a,  B)  = n 

o 

, - (w-X  ) /w 

kwe  p,  a.  < a < it 

p L — - 

(36) 

where  x = q (4)  - 

♦) 

(37) 

is  the  energy  gained  by  the  protons,  and 

a = sin 

The  limit  a.  is  given  by  (38)  when  a = a = 170°. 

L S p 

In  the  case  of  the  protons,  there  are  no  complications  due  to  reflec- 
tions because  both  the  electric  and  the  magnetic  forces  continually 
increase  the  proton  velocities  along  the  field. 


| (1  - J)  sin2  ag 
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~i  Ul 


w > X 


(38) 
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After  performing  the  pitch-angle  integration  in  (8) , the  number  density 
is  given  by  the  equation, 


P x/w 


III  no  6 P 
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where  g = — sin  a 
B p 

s r 


Note  that  for  B=B  , v = 0, 
s 


III/n  . n . 

n (B  ) = o(l  -cos  a ) 
P s -x-  p 


p ,a  as 


which  is,  as  expected,  equal  to  n (- — ),  the  number  density  due  to 

o 4it 

a Maxwellian  distribution  times  the  fractional  solid  angle  of  the 


source  cone.  By  equating  nI:[I(B  ) to  n1  + n11  = 1.711/cm3  at  B , the 
p p s 3 e e s 

value  of  n is  found  to  be  225.2/cm  . 
o 


2.3.4  Resultant  Potential  Distribution  and  Electric  Field 


The  number  density  of  the  protons  was  very  nearly  equal  to  the  total 

number  density  of  the  electrons  in  the  range  0.77  < B/B  < 1 for  the 

— s — 

potential  distribution  given  by  the  equation, 

♦,b)  - c.  dWBa~  °'77)2+  i]1'2  - n 

1 Ll 

where  = 6.494  x 10  ^ kV  and  = 2.999  x 10  \ By  differen- 
tiating (42)  with  respect  to  the  distance  along  the  magnetic  field, 
assuming  B to  vary  as  the  inverse-third-power  of  the  geocentric 
distance,  the  electric  field  is  found  to  be, 


p _ o/_L>2  B/Bg (B/Bo-0 . 77) 

^2  (Rs  + S)  «>(B)  + Cx) 


Here,  is  the  geocentric  distance  to  the  satellite  (13,660  km) 
and  S is  the  distance  above  the  satellite. 
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(39) 


(40) 


(41) 


(42) 


(43) 


I 

I 

The  potential  distribution  and  the  electric  field  are  plotted  in  Figure  8 

as  a function  of  B/B  . The  number  densities  of  the  protons,  n111,  and 
I SII  P 

the  electrons,  ne  + ne  » resulting  from  this  potential  distribution  are  , 

shown  as  a function  B/B  in  Figure  9.  Over  the  range  of  B/B  from  1 to 

s s | 

0.83,  the  number  densities  were  in  agreement  to  within  one  part  in  500.  j 

Toward  lower  magnetic  field  ratios,  the  number  densities  began  to  depart  j 

slowly  with  the  number  density  of  the  electrons  exceeding  that  of  the 
protons . 

As  mentioned  earlier,  the  end  points  of  the  potential  distribution  were 
fixed  at  those  values  given  by  the  electron  pitch-angle  distributions 
(0.8  kV  at  B^  and  0 kV  at  B^/B^O . 77)  . Only  the  shape  of  the  distribu- 
tion between  those  points  was  varied  in  this  analysis.  Nevertheless,  , 

since  the  number  densities  of  the  charges  are  quite  sensitive  to  the 
distribution,  it  seems  unlikely  that  appreciable  departures  from  those 

end-point  values  would  have  given  an  equally  good  agreement  of  the  , 

number  densities. 

■ 

The  measurements  of  Mozer  et  al.  (Ref.  3)  indicate  that  in  the  turbulent 
wave-field  region,  where  such  potentials  appear,  the  number  density  of 
the  thermal  plasma  may  be  as  high  as  50 /cm  with  fluctuations  of  about 
20%.  Data  on  the  plasma  density  for  the  event  reported  here  were  not 
obtained.  However,  even  if  the  thermal  plasma  density  was  high,  the 
preponderance  of  these  particles  as  discussed  in  Section  2.1  are  expected 
to  be  constrained  from  moving  along  the  magnetic  field  by  the  local  AC 
fields  of  the  waves,  which  are  tens  of  mV/m.  The  electric  field  shown 
in  Figure  8 is  less  than  about  1 mV/m.  Of  course,  the  runaway  electrons 
and  ions  do  move  rapidly  along  the  field,  but  the  majority  of  these 
particles  are  probably  included  in  the  distributions  used  in  this 
analysis . 

The  results  shown  in  Figure  8 together  with  those  in  Figure  6 indicate 
that  the  potential  distribution  has  a scale  comparable  to  an  earth's 
radius.  Furthermore,  the  measurements  of  the  Berkeley  group  (Ref.  16) 

on  this  event  of  the  turbulent  wave  field  and  the  field-aligned  current  j 

indicate  that  th?  potential  is  supported  by  anomalous  resistivity. 
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APPENDIX  2A 


CONDITION  FOR  REFLECTION  OF  ELECTRONS 
IN  ELECTRIC  AND  MAGNETIC  FIELDS 

If  the  energy,  w,  and  magnetic  moment,  y,  of  an  electron  are  known  at  B 
where  the  potential  is  q<j>=0,  the  parallel  component  of  the  electron  at 
h'  where  the  potential  is  ip'/q  is  given  by  the  equation, 


Wj^  = w - yB  + y(B-B')  - (ip-lp')  (44) 

Hwere,  w ” yB  equals  w^  at  B,  y(B-B')  is  the  increase  (or  decrease) 

in  w1  ^ due  to  the  change  in  B,  and  iJj-i p'  is  the  decrease  (or  increase) 

in  w^  due  to  the  change  in  the  potential. 


The  rate  of  change  of  w'  with  respect  to  B',  is,  therefore, 


V 


+ 


dip' 

dB 


(45) 


From  this  equation,  it  is  clear  that  for  upgoing  electrons  (B  decreas- 
ing), the  necessary  condition  for  mirroring  is  dw'/dB'  >0,  B'  > B', 

11  m 

where  B " is  the  mirror  point.  Thus,  from  (45) 
m 


> y , B ' > B_* 

QD  m 


(46) 


Similarly,  for  downgoing  electrons,  the  necessary  condition  for  mirror- 
ing is, 


3f  < - • »'  < 


(47) 


Since y = w sin  a/B,  the  pitch  angle  at  B for  mirroring  at  B',  following 

m 

from  Eq.  (44)  for  w.^=  0,  is  given  by  the  equation, 


.2  B . ip-<p  . 

sin  a = — ^ (1  - — ) 
c B w 

m 


(48) 
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The  extreme  value  of  a when  plotted  against  B'  occurs  at  the 

c m 

value  of  B'  given  by  the  equation, 

dil  , w (1  _ iztl  x = „ 


A second  differentiation  of  (46)  reveals  that  at  that  point,  a is  a 

c 

minimum  or  a maximum  depending  on  whether  the  curvature  of  \p', 

2 2 

d ii'/dB'  , is  positive  or  negative,  respectively, 
m 


As  discussed  in  the  text,  potential  functions  with  positive  curvature  were 

were  found  to  give  appropriate  results  for  the  B range  of  interest. 

Hence,  a was  a minimum  at  the  value  of  B'  given  by  the  solution 
c m 

of  (49). 
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Section  3 


RELATIONS  BETWEEN  GEOPHYSICAL  PHENOMENA  AND  IONOSPHERIC  SCINTILLATIONS 

3.1  INTRODUCTION 


The  WIDEBAND  satellite  program,  in  operation  now  for  over  a year,  has 
amassed  a large  body  of  data  on  the  occurrence  and  morphology  of  iono- 
spheric scintillations.  The  experiment  consists  mainly  of  a transmitter 
on  a near-polar  orbit,  broadcasting  on  ten  channels  from  VHF  to  S band 
[ref.  17].  The  signals  are  received  by  stationary  receivers  in  equatorial, 
mid-latitude,  and  auroral  regions.  The  major  shortcoming  of  the  program 
is  the  absence  of  energetic  particle  detectors  on  the  satellite.  It  has 
become  clear  from  coordinated  experiments  with  the  Chatanika  radar, 
rocket-borne  instruments,  and  ground-based  photometers  [ref.  18],  that 
there  is  a strong  association  between  certain  types  of  scintillation 
activity  and  particle  precipitation.  Efforts  are  underway,  too,  (by 
E.  J.  Fremouw  and  R.  R.  Vondrak  of  the  Stanford  Research  Institute,  and 
R.  D.  Sharp  of  Lockheed)  to  analyze  the  data  from  the  occasional  occur- 
rences where  the  WIDEBAND  data  could  be  correlated  with  particle  data 
from  Lockheed  instruments  on  the  S3-3  satellite.  In  this  report  we  will 
discuss  some  of  the  features  of  the  scintillations  and  their  apparent 
relationship  to  particle  precipitation.  A detailed  analysis  of  selected 
individual  cases  is  beyond  the  scope  of  the  present  work  and  will  no 
doubt  be  treated  by  the  above-cited  workers. 

3.2  MID-LATITUDE  OBSERVATIONS 

The  first  of  many  surprises  in  the  WIDEBAND  program  was  the  high  level  of 
scintillation  activity  at  mid-latitudes  [ref.  17].  The  receiver  station  at 
Stanford,  California,  was  originally  intended  for  checking  out  the  opera- 
tions; the  location  was  selected  solely  because  of  its  accessibility  to  the 
SRI  laboratories.  The  period  immediately  after  commencement  of  operations 
was  very  active  (see  Figure  20)  so  many  good  scintillation  observations 
were  observed  at  Stanford  during  the  first  two  months  [ref.  1].  However, 
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after  four  months  of  operation  the  Stanford  receiving  station  was  moved 
to  an  equatorial  location.  As  explained  in  our  last  DNA  report  [ref  l], 
not  enough  data  were  accumulated  at  Stanford  for  statistical  correlations 
with  geophysical  phenomena.  We  found  no  correlation  with  any  of  the 
important  geomagnetic  activity  indices. 

It  is  plausible  that  the  mid-latitude  scintillations  are  associated  with 
precipitation  from  the  trapped  radiation  belts,  but,  due  to  the  lack  of 
coordinated  particle  data,  it  has  not  been  possible  to  test  this  supposi- 
tion. The  sporadic  nature  of  the  mid-latitude  scintillations  could  be 
attributed  to  a number  of  possible  causes.  One  factor  that  determines  the 
spatial  structure  of  the  precipitation  regions  is  the  asymmetry  of  the 
geomagnetic  field;  drifting  electrons  tend  to  be  precipitated  in  distinct 
longitude  regions  (e.g.,  the  South  Atlantic  anomaly)  [ref.  1].  The  other 
important  factor  is  the  spatial  distribution  of  the  wave  sources  responsible 
for  dumping  electrons  from  the  radiation  belts.  It  has  been  suggested 
recently  that  man-made  VLF  emissions  are  very  effective  in  causing  energetic 
electron  precipitation  in  mid-latitudes  [ref  19].  Many  of  the  naturally 
generated  ELF/VLF  emissions  also  show  distinct  preferences  for  certain 
local  times.  It  may  be  concluded  that  more  scintillation  data  are  needed 
in  mid-latitudes  before  meaningful  statistical  associations  can  be  deduced. 

It  was  our  hope  that  high  latitude  wideband  receiver  stations  (Poker  Flat, 
Alaska)  could  be  used  to  back  up  the  mid-latitude  observations  by  observing 
scintillations  far  to  the  south,  in  the  mid-latitude  regions.  If  the 
propagation  path  is  assumed  to  intersect  the  scintillation  region  at 
300  km  altitude  in  the  F-region  ionosphere,  the  latitudinal  coverage  can 
readily  be  deduced.  The  limits  of  the  coverage  at  Poker  Flat  were  from 
about  50°N  geomagnetic  latitude  to  beyond  80°  geomagnetic  latitude.  However, 
the  scintillation  indices  increase  directly  with  the  total  electron  content 
along  the  propagation  path.  This,  together  with  the  difficulty  of  distinguish- 
ing spatial  structures  near  the  horizon,  limited  the  southward  coverage  to 
about  58°N  geomagnetic,  which  is  near  the  northern  limit  of  the  plasmasphere . 
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It  was  therefore  not  practical  to  use  the  Poker  Flat  data  for  routine  observa- 
tions of  mid-latitude  regions.  It  does  appear  that  the  mid-latitude  scintil- 
lations observed  at  Poker  Flat  were  similar  in  character  to  those  observed  at 
Stanford.  There  was  no  readily  discernible  pattern  that  might  be  correlated 
with  geomagnetic  activity. 

3 . 3 HIGH-LATITUDE  OBSERVATIONS 


The  scintillation  records  show  dramatic  changes  in  character  as  they  pass 
from  mid-latitude  through  the  plasma  trough  into  the  auroral  zone.  Figures 
10  through  19  [courtesy  C.  L.  Rino  of  Stanford  Research  Institute]  show 

some  sample  scintillation  records  from  the  Poker  Flat  station;  they  are  not 
necessarily  typical,  but  were  selected  to  illustrate  certain  features.  Figure 
10  shows  a quiet  pass.  Starting  at  the  right  there  is  perhaps  some  activity 
at  mid-latitudes  in  both  phase  (solid  line)  and  amplitude  scintillations 
(dashed  line) . There  is  only  a weak  correlation  (if  any)  between  the  phase 
and  amplitude  scintillations  on  the  right  sides  of  these  figures;  that  may 
be  true  of  most  of  the  mid-latitude  observations.  The  second  part  of  the 
figure  shows  the  relative-phase  scintillations  in  VHF,  UHF,  and  L band; 
this  record  is  of  value  mainly  in  indicating  where  the  VHF  records  are  sus- 
pect, as  at  the  right  side  of  Figure  10.  The  bottom  part  of  the  figure 
shows  the  tot-.'  electron  content  integrated  along  the  propagation  path.  The 
deep  minimum  shows  clearly  the  location  of  the  plasma  trough  which  is  expected 
at  the  southern  edge  of  the  auroral  zone.  The  record  in  Figure  10  shows 
weak  enhancement  of  the  electron  density  on  the  poleward  edge  of  the  trough, 
and  a slight  amount  of  scintillation  activity  probably  associated  with  aurorae 
beyond  63°  geomagnetic  latitude.  One  of  the  scintillation  spikes  here  is 
probably  spurious  because  the  propagation  path  is  nearly  aligned  with  the 
field  lines;  this  condition  is  usually  noted  by  a tick  mark  on  the  plot  of  BP 
angle,  here  at  17.7°.  At  low  and  moderate  levels  of  scintillation  activity 
the  activity  is  concentrated  near  the  poleward  edge  of  the  plasma  trough 
and  is  well  correlated  with  local  geomagnetic  indices  [communication  from 
C.  L.  Rino].  There  also  seems  to  be  some  correlation  between  the  phase 
and  amplitude  scintillation  indices  for  this  type  of  activity.  Figure  11 
shows  an  active  pass  and  illustrates  singular  features:  plasma  trough  with 
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Figure  10:  Scintillation  records  for  WIDEBAND  pass  7-2,  1105UT/8/6/76.  The  amplitude 
scintillations  are  denoted  S4^  (dashed  lines)  and  the  phase  scintillations  PRMS  (solid 
lines).  The  center  plot  shows  the  ratios  of  VHF/UHF/L-BAND.  The  lower  plot  is  of  total 
(integrated)  electron  content. 
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Figure  16.  Scintillation  records  for  WIDEBAND  pass  5-50,  0835UT/6/11/76. 


8 


activity  concentrated  at  its  poleward  edge.  The  phase  and  amplitude 
scintillation  indices  seem  to  follow  the  same  trend,  though  the  phase 
scintillations  outweigh  the  amplitude  scintillations  between  64°N  and  67°N 
geomagnetic  latitude.  For  convenience,  let  this  type  of  "typical"  scin- 
tillation activity  be  designated  Type  I. 

A morphological  survey  of  the  first  3-1/2  months  of  Poker  Flat  data 
reveals  that  there  may  be  yet  another  type  of  high  latitude  scintillation 
activity  (in  addition  to  Type  I and  the  characteristic  equatorial  and  mid- 
latitude activity).  These  events,  let  them  be  called  type  II,  are  charac- 
terized as  being  mainly  isolated  patches  or  bands  of  intense  phase  scintil- 
lations accompanied  by  weak  amplitude  scintillations.  They  usually  occur 
to  the  north  edge  of  the  auroral  zone  at  70°N  to  80°N  longitude.  Figures 
12  through  19  show  some  representative  cases  for  varying  levels  of  general 
activity.  Figure  12  is  a good  example  of  a mostly  quiet  pass  with  an 
isolated  phase  scintillation  spike  at  75°N.  (The  spike  at  64°  should  be 
ignored.)  There  is  no  corresponding  enhancement  of  the  amplitude  scintil- 
lations. The  total  electron  content  trace  shows  a patch  of  enhanced  ioniza- 
tion at  the  same  location  as  the  scintillation  spike,  indicating  that  there 
was  particle  precipitation  at  that  location.  The  electron  density  seemed 
to  be  elevated  over  the  entire  auroral  zone  from  75°N  to  64 °N.  This  might 
have  been  due  to  earlier  precipitation  that  had  ceased  by  the  time  of  the 
WIDEBAND  observation.  Unfortunately,  this  supposition  cannot  be  tested  be- 
cause of  the  absence  of  particle  precipitation  data.  Perhaps  the  chances 
of  correlating  WIDEBAND  observations  with  precipitation  (or  known  precipi- 
tation features)  will  improve  when  more  data  become  available  from  the 
winter  period,  1976-1977.  During  the  winter  it  is  feasible  to  observe  the 
aurorae  with  all-sky  cameras  and  photometers;  the  location  of  auroral 
arcs  and  precipitation  regions  can  be  readily  deduced  from  optical  observa- 
tions . 

Figures  13  through  15  show  additional  isolated  phase  scintillations  asso- 
ciated with  enhanced  F-region  electron  density  and  sometimes  with  weak 
amplitude  scintillations.  The  peak  at  the  far  left  on  Figure  14  is  an 
instrumental  effect  and  should  be  ignored. 
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Figure  16  is  distinguished  from  the  previous  examples  by  the  activity  over 
the  entire  latitude  range,  and  by  the  broad  band  of  elevated  total  elec- 
tron content.  The  plasma  trough  is  not  immediately  discernible,  but  there 
seems  to  be  a fall-off  in  TEC  far  to  the  right,  as  if  the  trough  had  been 
moved  to  L =3.  An  examination  of  the  geomagnetic  indices  for  0835  UT  June 
11,  1976,  when  the  satellite  was  overhead,  reveals  that  there  was  a major 
geomagnetic  storm  on  that  day  with  a sudden  commencement  at  0000  UT.  This 
was  one  of  the  largest  storms  during  the  entire  year;  and  by  0835  UT  the 
DST  value  had  reached  -51,  indicating  that  a fully-developed  ring  current 
had  been  established.  The  local  time  at  Poker  Flat  is  10  hours  earlier 
than  UT,  or  2235  LT.  The  WIDEBAND  observations  of  Figure  16  were  made 
just  before  local  midnight,  so  it  is  reasonable  to  suppose  that  the  scin- 
tillations at  79°N  and  72°N  were  related  to  intense  auroral  arcs.  The  ob- 
servations were  made  near  midsummer  so  the  auroral  activity  could  not  be 
observed  against  the  brightly  lit  sky. 

The  succeeding  WIDEBAND  passes  also  furnished  good  data,  shown  in  Figures 
17  and  18.  Figure  18  is  for  a pass  near  local  midnight;  overhead  at  0030 
LT.  This  time  was  coincident  with  the  maximum  D of  -61.  The  trough  had 
moved  to  its  southward  limit,  and  the  center  of  the  scintillation  activity 
had  moved  in  to  the  edge  of  the  trough.  The  very  strong  scintillations 
north  of  the  trough  are  not  of  the  expected  Type  I,  but  appear  to  be  almost 
entirely  phase  scintillations.  By  the  next  pass.  Figure  18,  at  0215  LT 
the  scintillation  activity  has  decreased,  though  the  phase  scintillation 
component  still  dominates.  The  TEC  curve  has  become  smoother,  indicating 
that  the  intense  precipitation  had  subsided.  But  the  tough,  and  the  Type 
II  scintillations  at  its  edge,  has  only  begun  to  move  back,  with  the  begin- 
ning of  the  recovery  phase  (DST=  ~52) . 

So  far,  efforts  to  correlate  the  Type  II  scintillation  with  S3-3  and  other 
satellite  data  have  been  unsuccessful.  There  were  several  storms  and  geo- 
magnetically  active  periods  shortly  after  the  S3-3  satellite  was  launched 
on  July  15,  1976  but  so  far  we  have  not  found  satisfactory  coordinated 
passes.  We  have  also  looked  at  ephemerides  of  other  satellites  in  the  first 
four  months  of  WIDEBAND  operations  and  found  no  suitable  coincidences. 
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3,4  MORPHOLOGY  OF  SCINTILLATIONS 


That  there  exist  several  distinct  types  of  high  latitude  scintillation 
activity  is  clear.  It  is  tempting  to  associate  the  Type  II  events  with 
auroral  arcs  and  precipitation  regions.  The  cases  described  are  certainly 
in  the  expected  locations.  The  associated  TEC  enhancements  are  suggestive 
but  not  conclusive.  In  the  preceding  figures  the  bumps  in  the  TEC  were 
small,  and  did  not  indicate  large  increases  in  the  electron  densities. 
Figure  19  perhaps  illustrates  the  electron  density  increases  more  clearly. 
This  remarkable  pass  also  coincided  with  the  main  phase  of  a geomagnetic 
storm  (0135  LT,  Sept.  2,  1976,  Dg^  = -31  ).  The  TEC  excursions  here 
indicate  order-of-magnitude  increases  in  the  electron  density,  again  asso- 
ciated with  high  phase  scintillation  activity. 


A more  complete  picture  of  the  morphology  of  the  scintillations  is  pre- 
sented in  Figure  20.  Here  only  the  widths  of  the  rms  phase  scintillation 
peaks  that  exceeded  1 radian  were  plotted.  Three  degrees  of  shading  are 
used  to  indicate  levels  of  activity  (these  levels  are  not  the  same  as  the 
activity  levels  described  in  reference  17) . No  attempt  was  made  to  separate 
out  the  Type  I activity  from  the  Type  II  activity.  Also  shown  is  the  loca- 
tion of  the  plasma  trough,  estimated  from  the  TEC  records.  The  pass  num- 
bers and  overhead  times  are  listed  for  all  the  WIDEBAND  nighttime  passes 
up  to  Sept  3,  1976.  To  the  far  right  are  listed  several  global  geomagnetic 
indices.  The  daily  sums  of  K are  weighted  toward  the  polar  regions,  and 
indicate  the  level  of  storm  and  substorm  activity.  The  hourly  D records 
are  weighted  toward  the  equator,  and  are  valuable  for  assessing  the  develop- 
ment of  the  ring  current  (at  L -3-5)  during  magnetic  storms. 


Several  features  stand  out  in  Figure  20.  The  centers  of  scintillation  ac- 
tivity and  the  plasma  trough  move  northward  and  southward  in  response  to 
geomagnetic  activity.  The  storms  of  June  11,  June  30,  and  Sept.  2 stand  out 
especially.  But  the  most  striking  aspect  of  the  morphology  is  the  regularity 
of  the  structures  of  roughly  1-2°  latitudinal  extent.  Translated  into  hori- 
zontal distance  these  structural  features  are  of  the  order  of  100  to  200 
km  wide.  The  structure  is  real;  the  spatial  sampling  rate  is  much  smaller 
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Figure  20a.  The  morphology  of  high-latitude  WIDEBAND  scintillation  observations  up 
to  9/3/76.  The  shaded  bars  are  regions  of  intense  phase  scintillations . The  location 
of  the  plasma  trough  (heavy  solid  line)  is  estimated  from  the  total  electron  content 
as  deduced  from  the  WIDEBAND  data. 
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Figure  20b.  The  morphology  of  high-iatitude  WIDEBAND  scintillation 
observations  up  to  9/3/76.  (continued) 
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than  the  smallest  structural  features  indicated  on  the  figure.  The 
structural  patterns  are  consistent  with  their  interpretation  as  auroral 
features.  When  the  scintillation  peaks  are  plotted,  as  here,  the  iso- 
lated peaks  stand  out  quite  clearly.  Taken  as  a whole,  the  morphology 
of  the  scintillation  regions  very  much  resembles  the  morphology  of  precipi- 
tation events  and  aurorae  in  a meridian  near  local  midnight. 

3.5  DISCUSSION  OF  RESULTS 

Our  tentative  interpretation  of  the  Type  II  scintillations  is  that  they  are 
related  to  the  large  scale  electric  field  structures  discussed  in  the  pre- 
ceding sections.  Their  statistical  occurrence  is  similar  enough  that  more 
work  should  be  done  to  find  and  interpret  coordinated  observations.  The 
bodies  of  WIDEBAND  and  S3-3  data  are  large  enough  that  good  coordinated 
pass  should  be  found  in  the  later  records  that  haven't  yet  been  analyzed. 

It  must  be  recognized,  however,  that  a triple  coordination  is  needed: 
between  two  satellites  and  an  active  region.  Ideally,  these  questions  could 
best  be  answered  if  a successor  to  the  WIDEBAND  satellite  were  put  in  orbit 
with  on-board  particle  detectors  that  could  delineate  the  precipitation 
features  and  the  electric  field  structures. 

Why  the  Type  II  scintillations  should  be  associated  with  electron  precipi- 
tation is  not  well  understood.  The  scintillations  are  due  mainly  to 
processes  in  the  F region,  at  altitudes  near  300  km.  The  normal  type  of 
scintillation,  presumably  Type  I,  is  expected  to  result  from  small  scale 
field  aligned  irregularities  formed  by  plasma  instabilities.  The  classical 
example  is  in  the  equatorial  electrojet,  where  streaming  type  instabilities 
cause  the  formation  of  field-aligned  tubes  with  enhanced  densities.  The 
characteristics  of  the  regions  responsible  for  type  II  scintillations  may 
become  better  known  when  some  cases  are  examined  in  detail,  to  ascertain 
the  statistics  of  the  scintillation  components  [see  ref.  17].  At  the 
moment  the  only  solid  pieces  of  evidence  are  their  size  and  occurrence, 
and  the  predominance  of  phase  scintillations.  All  these  characteristics 
tend  to  contradict  the  interpretation  of  Type  II  as  being  due  to  field- 
aligned  plasma  irregularities.  One  possible  interpretation  is  that  the 
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Type  II  scintillations  result  when  the  propagation  path  passes  through  the 
"pencils"  of  ionization  caused  by  the  precipitating  electrons  that  are 
responsible  for  the  fine  scale  auroral  features.  Optical  observations 
have  not  set  lower  limits  on  the  sizes  of  auroral  forms,  they  might  have 
sizes  down  to  the  gyroradii  of  the  electrons.  A 12  keV  electron,  for 
instance,  has  a gyroradius  in  the  F-region  of  about  1 meter.  This  gives 
a scale  size  rather  larger  than  the  anticipated  irregularities;  somewhat 
smaller-scale  irregularities  with  a fairly  high  degree  of  spatial  organi- 
zation would  be  needed  to  cause  intense  amplitude  scintillations,  especially 
in  the  UHF  and  higher  frequencies. 
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Section  4 

PLASMA  IRREGULARITIES  RESULTING  FROM  A H.A.N.E. 


4.1  INTRODUCTION 


( 
1 

On  the  other  hand,  the  hot  plasma,  i.e.,  the  nuclear  debris  and  the  relati-  ! 

vistic  electrons,  injected  into  the  magnetosphere  produces  electric  fields 
and  currents  which  extend  to  thousands  of  kilometers  both  along  and  across 
magnetic  field  lines.  These  currents  are  susceptible  to  streaming  insta- 
bilities which  produce  electrostatic  waves  and  associated  irregularities. 

The  field-aligned  currents  can  produce  irregularities  to  altitudes  of 
several  thousand  kilometers,  as  discussed  in  Section  2.1  (see  Ref.  7); 
the  cross-field  currents  produce  the  irregularities  in  the  ionosphere, 
but  they  also  propagate  to  high  altitudes  along  the  magnetic  field  (Ref. 

20  and  21).  All  these  irregularities  resulting  from  the  hot  plasma  extend 
to  much  greater  distances  than  those  due  to  the  locally-heated  atmosphere, 
but  their  duration  may  not  be  as  long. 

As  discussed  below,  the  injected-debris  plasma  produces  both  field-aligned 
and  cross-field  currents  that  may  be  unstable  to  the  growth  of  electro- 
static waves.  The  resulting  irregularities  are  expected  to  occur  north- 
ward of  the  burst  to  latitudes  near  the  auroral  zone.  The  relativistic 
electrons  emitted  by  the  fission  debris  are  expected  to  produce  only  cross- 
field currents  which  may  be  unstable.  However,  owing  to  the  eastward- 
drift  motion  of  these  electrons,  the  resulting  irregularities  are  expected 
to  occur  over  a considerable  distance  east  of  the  burst  meridian. 


A high-altitude  nuclear  explosion  deposits  considerable  energy  in  the  upper 
atmosphere  and  injects  hot  plasma  into  the  magnetosphere.  The  configura- 
tion of  the  thermal  plasma  due  to  the  enhanced  ionization  and  heating  of 
the  atmosphere  is  unstable.  It  expands  along  the  local  magnetic  field 
and  breaks  up  into  irregularities  which  have  a broadband  spatial  structure. 
The  gradient-drift  instability  appears  to  be  the  most  important  mechanism 
responsible  for  these  irregularities,  and  much  attentic  as  been  given 
to  the  assessment  of  effects  due  to  this  mechanism. 
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4.2  CURRENTS  DUE  TO  INJECTED  DEBRIS 


A large  fraction  of  the  debris  from  a H.A.N.E.  moves  upward  along  the  mag- 
netic tube  that  is  referred  to  as  the  debris  tube.  The  energy  of  the 
ionized  debris  is  high,  and  it  therefore  drifts  rapidly  toward  the  west 
owing  to  the  configuration  of  the  earth's  magnetic  field.  On  the  other 
hand,  the  energetic  electrons  emitted  by  the  debris  drift  toward  the 
east.  The  resulting  charge  separation  produces  a large  electric  field 
across  the  tube,  directed  toward  the  east.  This  field  causes  the  plasma 
in  the  tube  to  drift  outward  (E  X B direction)  toward  higher  L values. 
Simultaneously , rhermal  plasma  currents  flow  along  the  magnetic  field  lines 
containing  the  excess  charges  and  across  magnetic  field  lines  in  the 
ionosphere  tending  to  neutralize  the  excess  charges. 

In  an  investigation  of  the  debris  motion  (Ref.  22)  it  was  found  that  for 
bursts  in  the  megaton  range  the  excess  charges  are  never  fully  neutralized. 
Hence,  the  debris  tube  continues  to  drift  outward  to  the  limit  of  the 
trapping  region  of  the  magnetosphere.  Moreover,  the  investigation  reveals 
that  the  field-aligned  currents  exceeded  the  instability  threshold  value 
which  was  given  by  Kindel  and  Kennel  (Ref.  7)  and  verified  by  the  S3-3 
satellite  data  (Ref.  3).  Plasma  irregularities  similar  to  those  observed 
with  the  S3-3  satellite  - extending  to  altitudes  in  excess  of  about  8000 
km  - should  therefore  be  expected  along  all  "closed"  field  lines  north 
of  the  burst. 

The  electric  field  across  the  tube  also  appears  in  the  E-region  of  the 
ionosphere.  There,  the  electrons  drift  essentially  according  to  the  E X B 
drift  velocity,  but  the  ions  are  greatly  restricted  by  collisions  with 
the  neutral  constituents  of  the  atmosphere.  Hence,  the  relative  drift 
velocity  of  the  electrons  with  respect  to  the  ions  depends  on  the  magnitude 
of  the  electric  field.  According  to  the  Farley-Buneman  criterion,  the 
resulting  ionospheric  current  becomes  unstable  if  the  relative  electron- 
ion  drift  velocity  exceeds  the  thermal  velocity  of  the  ions. 
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Using  the  model  for  the  debris-tube  motion  described  in  Ref.  (22)  , the 
electric  field  in  the  ionosphere  was  computed  for  the  parameters  of  Star- 
fish and  the  USSR  burst  of  Oct.  22,  1962.  The  thermal  velocities  of  the 
ions  were  computed  assuming  that  they  were  in  thermal  equilibrium  with 
the  neutral  atmosphere  in  the  E-region.  It  was  found  that  the  E X B 
drift  velocity  of  the  electrons  exceeded  10^  m/sec  at  latitudes  higher 
than  19°N  in  the  case  of  Starfish.  For  the  USSR  burst,  the  drift  velocity 
exceeded  10~*  m/sec  at  latitudes  above  41°N.  The  relative  drift  velocity, 
therefore,  far  exceeded  the  ion  thermal  velocities,  which  were  found  to 
be  in  the  range  420-480  m/sec,  depending  on  diurnal  and  solar-activity  con- 
ditions. Hence,  as  stated  above,  plasma  irregularities  following  a high- 
yield  burst  should  be  expected  from  the  ionosphere  to  altitudes  of  8000 
km  or  more  on  field  lines  extending  from  the  burst  point  to  the  auroral 
zone. 


4.3  CURRENTS  DUE  TO  INJECTED  ELECTRONS 

As  the  debris  tube  drifts  outward,  the  electrons  emitted  by  the  debris 
continually  drift  out  of  the  tube  toward  the  east.  Owing  to  velocity 
dispersion  the  electron  distribution  spreads  toward  the  east  and  large- 
scale  electric  fields  ensue  from  the  excess  charges.  Again,  thermal- 
plasma  currents  along  magnetic  field  lines  tend  to  neutralize  the  excess 
charges  but  these  charges  cannot  be  neutralized  in  times  less  than  the 
time  required  for  a hydromagnetic  wave  to  traverse  a field  line.  During 
that  time  the  distribution  of  the  injected  electrons  changes  as  the  elec- 
trons drift  across  magnetic  field  lines  toward  the  east;  thus,  the  elec- 
tric field  persists  while  the  electrons  remain  bunched  during  their  drift 
motion. 

In  a model  discussed  in  Ref.  (22),  the  electric  field  across  the  magnetic 
field  due  to  the  injected  electrons  was  computed,  taking  into  account 
the  effects  of  the  changing  injected-electron  distribution  and  the  neutral- 
izing currents.  The  field  was  computed  for  the  electrons  injected  by  Star- 
fish and  the  USSR  test  of  Oct.  22,  1962.  Again,  this  field  was  transformed 


to  the  E-region  of  the  ionosphere,  and  the  resulting  E X B drift  velocity 
of  the  electrons  was  determined.  It  was  found  that,  for  the  Starfish 
distribution,  the  relative  drift  velocity  of  the  electrons  with  respect 
to  the  ions  exceeded  the  thermal  velocity  of  the  ions  in  the  latitude 
range  20-56°N  and  in  a longitudinal  interval  extending  from  the  burst 
meridian  to  about  60°E.  For  the  USSR  burst,  the  relative  drift  velocity 
exceeded  the  ion  thermal  velocity  at  latitudes  between  41°N  and  51°N, 
and  at  longitudes  extending  to  about  40°  east  of  the  burst  meridian. 

The  field-aligned  currents  due  to  the  injected  electrons  did  not  exceed 
the  Kindel  and  Kennel  instability  criterion.  Hence,  only  the  ionospheric 
irregularities  should  be  expected  over  the  regions  given  above. 

The  duration  of  the  irregularities  is  uncertain.  It  is  likely  that  they 
will  persist  at  least  as  long  as  the  driving  perturbation  remains.  For 
the  case  of  the  debris  motion,  che  perturbation  may  last  only  a few  hun- 
dred seconds.  For  the  case  of  the  injected  electrons,  it  lasts  about  1/2 
hour  at  L : 2.  However,  the  duration  of  the  electric  field  is  proportional 
to  the  longitudinal  drift  period  of  the  electrons.  Hence,  since  the 
drift  period  is  inversely  proportional  to  L,  the  duration  of  the  pertur- 
bation should  decrease  toward  higher  L values. 
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